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SUMMARY 


A system of partial differential equations of motion for a steadily rotat- 
ing circular ring is derived- The ring is homogeneous with cross-sectional 
dimensions which are small relative to the overall diameter. Small deflections 
from the steady- state condition are assumed so that the equations of motion are 
linear. 

Four sets of structural modes are examined. The torsional modes consist 
of rotations about the line of centroids of the cross section with small deflec- 
tions perpendicular to the plane of the ring. The out-of-plane bending modes 
consist of deflections perpendicular to the plane of the ring accompanied by 
rotations about the line of centroids. The in-plane bending modes consist of 
radial deflections with small tangential motions. The compression modes con- 
sist of stretching and compression of the line of centroids with small radial 
deflections . 

For each mode, the motions are written as sinusoidal functions of time and 
position around the circumference of the ring. Several low- numbered modes are 
shown to be rigid body motions which agree with commonly known results of rigid 
body mechanics- Formulas for the frequencies of vibration and the ratios 
between different perturbation variables are derived. In many cases the for- 
mulas are simplified- For example, for low spin rates the frequencies of vibra- 
tion of the bending modes approach the frequencies for a linear beam. At high 
spin rates, these frequencies approach the frequencies for a string under 
tension . 

Several figures are used to demonstrate the spatial and temporal charac- 
teristics of the structural modes. The derived formulas are applied to an 
example problem. The results of the example are plotted to illustrate the 
effects of spin rate and mode number upon the frequencies of vibration and 
other parameters which quantify the vibrational characteristics of the spinning 
ring . 


INTRODUCTION 

Attitude control of spacecraft is accomplished by generating torques, 
either externally as may be done by small rocket thrusters or internally by 
stored momentum control devices. One example of a stored momentum control 
device is the control moment gyro (CMG) . (See ref. 1.) The CMC consists of a 
momentum wheel which is supported in a gimbal system. Torques applied to the 
support bearings of the wheel cause it to process. The precession is used to 
change the spacecraft attitude. In reference 1, three CMG*s are arranged so 
that the nominal spin axis of each CMG is aligned with one of the spacecraft 
axes. Another example of a stored momentum control device is the so-called 
scissors gyro concept described in references 2 and 3. In this approach, two 


momentum wheels are arranged so that their nominal spin axes coincide but their 
spin directions are opposite. Precession of these wheels in opposite directions 
results in a net transfer of angular momentum to the spacecraft. 

The effectiveness of a momentum wheel in a given application depends on 
the amount of angular momentum which can be produced about its spin axis. The 
moment of inertia of the wheel and its operating spin rate determine the angu- 
lar momentum. In general, devices which are to be placed in spacecraft are 
designed with the least possible mass so that required launch system capabili- 
ties are minimized. Meeting the inertia requirement of a momentum wheel while 
minimizing its mass leads to the annular momentum control device (AMCD) . (See 
ref- 4.) The AMCD consists of a ring (annulus) which is supported by noncon- 
tacting magnetic bearings and is spun by a noncontacting electromagnetic motor. 
For large precessions, the ring and its magnetic bearing assembly can be 
rotated in a gimbal system. Small accurate precessions may be obtained by 
varying the magnetic forces in the bearings. An advantage of the AMCD concept 
over conventional shaft-mounted momentum wheels is that the mass of the center 
portion of the wheel is eliminated with a relatively small reduction in the 
moment of inertia. 

With the possibility of future spacecraft being gigantic, AMCD*s of very 
large dimensions can be envisioned. As AMCD*s increase in overall diameter and 
decrease in cross-sectional area, their flexibility increases. Thus, structural 
modes in the device must be considered in the design . 

Several analyses of structural vibrations occurring in rotating systems are 
in the literature- For example, reference 5 examines the flexural vibrations of 
rotating shafts with various loads, support methods, and unbalance conditions. 
This work treats high-speed rotary systems in which the shafts rotate at speeds 
in excess of their critical speeds. When the rotational period approaches the 
period of transverse flexure modes of the shaft, divergent oscillations occur 
and cause severe bearing vibrations and possible structural failure. A second 
type of instability occurring in rotating systems is discussed in reference 6. 
This work analyzes a self-excited oscillatory helicopter rotor instability. 

This potentially dangerous phenomenon was shown to be caused by a conversion of 
the rotational energy of the rotor into oscillatory energy of the blades. Both 
of the above types of instability can be present in centrifuges. Reference 7 
examines the stability characteristics of a centrifuge which is mounted in a 
large spacecraft. The shaft cri tical-speed instability arises because the mass 
of the centrifuge is concentrated along the arm connecting the gondala and 
counterbalance with the hub. The second oscillatory instability is associated 
with flexure of the arm and movement of the hub. The speed ranges within which 
these instabilities occur can be adjusted by altering the mass distribution, 
stiffness properties, and damping characteris tics of the centrifuge. 

Circular rings and arches have long been employed in construction (build- 
ings, aqueducts, bridges, etc.) and in cylindrical vessel designs (tanks, pipe- 
lines, aircraft, etc.). The strength, rigidity, and possibility of buckling 
are considerations in these designs. Calculated static stresses, deflections, 
and elastic stability characteristics of circular rings and arches with various 
loading conditions can be found in reference 8. A brief analysis of the 


stresses within a steadily rotating hollow disk can be found in reference 9. 

This work does not include vibrations or structural modes. 

An early analysis of the vibrations of a circular ring can be found in 
reference 10. This analysis considers a slender ring which has a circular 
cross section, is not spinning, and has no external influences (forces) acting 
upon it. The linear analysis produces formulas for the frequencies of oscilla- 
tion for flexural vibrations in the plane of the ring, flexural vibrations per- 
pendicular to the plane of the ring, and extensional vibrations along the line 
of centroids of the ring. Reference 11 examines the forced response character- 
istics of the flexural vibrations perpendicular to the plane of the ring. 

Forcing terms are added to the analysis of reference 10 and generalized coordi- 
nates are used in developing the response of the vibrational modes to externally 
applied forces. A corresponding analysis of the forced response characteristics 
of the flexural vibrations and the extensional vibrations in the plane of the 
ring is also developed in reference 11. The analysis of the flexural vibrations 
is validated by a physical experiment. A nonlinear analysis of the flexural 
vibrations in the plane of the ring is developed in reference 12. Nonlinear 
resonance curves which demonstrate a nonlinearity of the softening type are 
developed. The analysis is supported by physical experiment. An analysis of 
the forced vibrations of a ring which is spinning is presented in reference 13. 
The linearized structural mode dynamics are developed in terms of generalized 
coordinates by means of Lagrange's equations. This work considers a thin, 
flexible ring with flexural and extensional vibrations in the plane of the ring 
and flexural vibrations perpendicular to the plane of the ring. 

The present paper derives linearized equations of motion of a spinning 
circular ring by application of Newton's law and Euler's equation following the 
method of reference 14. The resulting set of linear partial differential equa- 
tions describe the small- amplitude dynamics of the flexural and extensional 
motions in the plane of the ring, the flexural motions perpendicular to the 
plane of the ring, and the twisting motions about the centroid of the ring. 

The vibration frequencies and modal characteristics of the freely vibrating 
ring are examined in detail with emphasis on physical interpretation of the 
mathematical results . 

The equations of motion are developed in appendix A. These equations are 
examined in the main text. The individual vibration modes are identified and 
simple expressions for their characteristic properties (for example, vibration 
frequencies) are derived. The spatial and temporal characteristics of each 
vibration mode are illustrated in several figures. Formulas developed in this 
paper are applied to an example ring and numerical results are presented. 


SYMBOLS 

2 

A area of ring cross section, m 

Aj , B j ,C j , . . . ,H j modal amplitudes of homogeneous solutions of wave equations 

ratio between amplitudes of twist angle and out-of-plane deflection, 
rad/m 


3 



/\ XS 

/\ /N ✓S, 

®r'®0'®z 


G 


Ir'Iz 

/\ /N 

i/ j /k 

J 

j 

M 

R 

AR 

R 


modulus of elasticity. Pa 

nondimens ional constant equal to EIj^/GJ 

unit direction vectors of curvilinear coordinate system 
unit direction vectors of cylindrical coordinate system 
force vector acting on cross section 

components of F in curvilinear coordinate system 
shear modulus. Pa 

ratio between amplitudes of tangential perturbation and radial pertur- 
bation, nondimensional 

polar moment of inertia of ring cross section, m^ 
area moments of inertia of ring cross section, m^ 
unit direction vectors of inertial coordinate system 
torsional constant of ring cross section, 
mode number 

moment vector acting on cross section 

components of M in curvilinear coordinate system 
radius of unstressed ring, m 
radius change caused by centrifugal force 
radius change caused by a small increment in spin rate 
location of ring centroid in inertial space 
location of ring fiber in inertial space 


nondimensional constant equal to \jj/AR^ 

s/ip/AR 


nondimensional constant equal to 


nondimensional constant equal to 
nondimensional constant equal to sJi^/ar^ 

arc length from reference point on ring to any other point on ring 
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# • - • 


polynomial coefficients 


t 

u 

u 

V 


W 

Wf f ^2 ' ^3 


X,Y,Z 


x,y, z 


x,y,z 


^f 


^in 

3 

Y 

£ 

^max 




9p'<i>p 


time, sec 

radial perturbation variable, m 
in-plane nodal frequency, rad/sec 
tangential perturbation variable, m 
in-plane nodal velocity, m/sec 

velocity of transverse sound waves, equal to yG/p, m/sec 
velocity of longitudinal sound waves, equal to '/e/P, m/sec 
externally applied control force density, N/m 

components of W in curvilinear coordinate system 
out-of-plane deflection perturbation variable, m 
inertial coordinate axes 

components of displacement for rigid body translation modes 
components of velocity for rigid body translation modes 
radial location of fiber with respect to centroid 
radial location of the innermost fiber with respect to centroid 
axial location of fiber with respect to centroid 
angle-of-twist perturbation variable, rad 
shearing strain, nondimensional 
strain, nondimensional 

maximum strain occurring among fibers caused by centrifugal force 
strain of line of centroids caused by centrifugal force 
angular components for rigid body rotation modes, rad 
angular rate components for rigid body rotation modes, rad/sec 
angular components for rigid body precession modes, rad 


X 


tangential distance from X-axis to fixed observation point 



p 


mass density of unstressed ring, kg/m^ 

0 stress. Pa 

0^ maximum allowable stress. Pa 

^max maximum stress occurring among fibers caused by centrifugal force, Pa 

T shearing stress. Pa 

Q. spin rate of ring, rad/sec 

spin rate required to produce the maximiom allowable stress CT^, rad/sec 

limiting spin rate for static stability, rad/sec 

0) frequency of oscillation, rad/sec 

velocity of transverse sound waves divided by the radius of the ring, 
rad/sec 

velocity of longitudinal sound waves divided by the radius of the 
ring, rad/sec 

Subscripts: 

B in-plane bending mode 

C in-plane compression mode 

D out-of-plane bending mode 

j mode number 

P progressive wave 

R regressive wave 

S relative to inertial space 

T torsional mode 

u in-plane mode 

w out-of-plane mode 
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Special symbols: 


A dot over a variable (e.g., v) denotes differentiation with respect to 

time . 


A prime after a variable (e.g., v*) denotes a spatial derivative with 

respect to s. 

Vertical bars enclosing a vector quantity (e.g., I I ) <^enote the length 
of the vector. 


EQUATIONS OF MOTION 

The linearized equations of motion for a steadily rotating slender ring 
are derived in appendix A. The properties of the ring and the assumptions used 
in the derivation are as follows: 

1. The ring is circular and has a uniform cross section. 

2. One principal axis of the cross section is parallel to the spin axis. 
The other principal axis lies in the plane of the ring. 

3. The ring is rotating about its axis at a constant rate. 

4. Perturbations of the ring from its steady state are small so that pro- 
ducts of the perturbation variables with each other may be neglected. 

5. Externally applied forces are small so that products of these forces 
with the perturbation variables may be neglected. 

6. The stress-strain relationships of the material of the ring are linear. 

7. The damping of the material of the ring is negligible. 

8. Deformations caused by transverse shear are ignored. 

9- Rotary inertia terms are included. 

10. The relative change in the radius of the ring caused by its steady 
rotation is small, 

11. The rate of rotation is small relative to the speed of sound traveling 
around the circumference of the ring. 

The small-perturbation equations of motion derived in appendix A 
(eqs. (A54) to (A57) ) are 
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( 1 ) 


u - 2f2v - - ojg^S^^ (R^u"" + 2R^u" + u) - ojg^(Rv' + u) 

+ f^^(R^u" - Rv') + S 2 ^(R^ii" - Rv') 


^2 0 

V + 2f2u = + (jOo^(Ru* 

pA ^ 

+ R^v" ) 

+ + Sg^CRvi' - V) 

(2) 

W = (rV’" 

pA s r 

- r33") 

+ (R^w" + R^3") 



- 2^3' 

- fiV') 

(3) 

Sp2 (3 + fiw' ) + _ 

S22)^(w' 

- f23) = 0)j.2sj2 (r 23" + Rw") 



+ a)g^Sr^(Rw" - 6) (4) 


where 

Q is spin rate. 

R is the radius of the unstressed ring. 

A is the area of the ring cross section, 
p is the mass density of the unstressed ring. 

u is a small-perturbation displacement in the radial direction. 

V is a small-perturbation displacement in the tangential direction. 

w is a small-perturbation displacement perpendicular to the plane of the 
ring. 

B is the small-perturbation angle of rotation about the line of centroids. 

, W 2 f and are external distributed force densities (force per unit 

length) acting in the radial direction, tangential direction, and per- 
pendicular to the plane of the ring, respectively. 

u, V, w, W 2 # and are functions of s and t. 

t is time. 

s is the arc length from a reference point on the ring to any other point 
on the ring. 
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A dot over a variable denotes differentiation with respect to time. 


A prime after a variable denotes a spatial derivative with respect to s. 
The constants appearing in equations (1) through (4) are defined as follows: 




0)v 



(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


where 

E is the modulus of elasticity. 

G is the shear modulus . 

Ip is the polar moment of the ring cross section. 

J is the torsional constant of the cross section. 

I^ and are the area moments of inertia of the cross section with 

respect to the radial direction and the spin axis, respectively. 

0)g is the speed of propagation of longitudinal sound waves divided by the 
radius of the ring. 

(jL)j. is the speed of propagation of transverse sound waves divided by the 
radius of the ring. 

Sj, Sp, S^, and are analogous to the inverse of the slenderness 

ratio used in the design of columns . 


More specific information concerning the definition of the variables and con- 
stants used above may be found in appendix A. 
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STATIC STABILITY 

This section examines the possibility of a static instability occurring in 
the spinning ring. Given an increase in the radius of the ring Ar, if the 
tension change required to balance the centrifugal force is greater than the 
increase in tension caused by stretching, then the radius will grow without 
bound. Given the unperturbed spinning ring, the tension caused by stretching 
derived in appendix A is 


F 


2 


EA 1 + 


\ AR 

ar2;r 


(11) 


Also from appendix A (eq. (A48) ) , the tension required to balance centrifugal 
force is 


F 


2 


pA 1 



( 12 ) 


The equilibrium strain in the ring is obtained by equating equations (11) 
and (12): 




AR 

R 


(Rfi) ' 




AR^ 


{bU} ■ 


(13) 


The limiting spin rate below which the ring is statically stable and above 

which the ring is unstable, is obtained by differentiating equations (11) 
and (12) with respect to the change in radius and equating the results: 





(14) 


Therefore, the instability occurs when the spin rate exceeds the speed of sound 
in the material of the ring divided by the radius of the unstressed ring. 


To determine whether this limiting spin rate would be achieved in practice, 
the spin rate required to produce the maximum allowable stress in the 

material can be computed and compared with The maximum stress in the 

steadily spinning ring is obtained from equations (All) and (A12) of appendix A: 


^max ^^max 


E£ 


o 


1 + 


^in 

R 


(15) 
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mil II I nil I 


mil 



where the perturbation variables are all zero and ^min radial distance 

between the innermost fiber and the centroid and has a negative value. Substi- 
tution of equation (13) into equation (15) , letting the maximum strain equal the 
allowable strain, and solving for the spin rate yield the maximum allowable spin 
rate based on the maximum allowable strain and the dimensions of the ring: 





1 + 


AR" 


1 + 


^min 

R 




1 + 


(16) 


For example, let the ring have a square cross section of dimension a on 
each side with a < R/10. For this case. 


a2 1 

^ ^ i — << 1 


AR^ 12R^ 


(17) 


^2^ = « 1 
R 2R 20 


(18) 


Also, consider the ring to be fabricated from steel for which (ref. 15) 


~ 0.7 X 10 ^ « 1 


(19) 


Applying the approximations of equations (17) , (18) , and (19) to equation (16) 

yields an approximate formula for the allowable spin rate based on the allowable 
stress of the material: 




a 



( 20 ) 


Comparing equations (20) and (14) , 



(21) 


which implies that the ring would experience structural failure well before the 
limit of equation (14) . This result applies to all relatively rigid materials 
for which E >> O^. 


11 



OUT-OF-PLANE HOMOGENEOUS SOLUTIONS 


The out-of-plane vibration inodes (deflection perpendicular to the plane of 
the ring and twisting about its centroid) are examined in this section. The 
equations of motion of the out-of-plane modes are given by equations (3) 
and (4) . In the absence of external influences (control forces and gravity ) , 
the distributed force acting perpendicular to the plane of the ring W 3 is 
zero. Because of the geometry of the ring/ the structural perturbation vari- 
ables are periodic in the variable s, that is. 


w (s 

+ 2ttR) = w(s) 

( 22 ) 

6 (s 

+ 2frR) = 6 (s) 

(23) 


for all s. As in all wave equations, there is a family of solutions to equa- 
tions (3) and (4) which satisfy the constraints of equations (22) and (23) . 
Assume that each mode has the solution 


w(s, t) 
3(S/t) 


Aj cos 


C . cos 
1 


if * “wj') * (¥ * “Wj‘) 

{f * “wj") * “j “»:•=) 


(24) 


(25) 


where the mode number j is a nonnegative integer and 03^j is the vibration 
frequency. 


Substitution of equations (24) and (25) into equations (3) and (4) with 
= 0 yields the following algebraic equations : 


(1 + _ j2(j2E^ + - (1 - 

= [j^R(EQ + l)Wr^Sj2 - (26) 
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where 


Eg = 


EI, 


GJ 



The constant Eq relates the out-of-plane bending stiffness to the torsional 
stiffness, and is the ratio between the amplitudes of the angle of rota- 

tion of the cross section and the deflection of the centroid perpendicular to 
the plane of the ring. 

Solving equations (26) and (27) by eliminating yields a polynomial 

in the out-of-plane vibration frequency : 

= 0 (28) 


where 


T 


a 


T 


b 


[(i - * j2(i . 

+ - Sj-^ + j^(sp^ + 2Sp4jjs^2 


= 4j^(l + 



The roots of equation (28) are the vibration frequencies of the out-of-plane 
modes. The proportionality constants C^j between the twist and the out-of- 
plane deflection for each of these modes are obtained by substituting the values 
of the frequencies into either equation (26) or (27) . 

The remainder of this section first examines the vibration characteristics 
when the spin rate is zero. Then the modes with mode numbers of zero and one 
with nonzero spin rate are examined. Finally, the vibration characteristics of 
the higher numbered modes with nonzero spin rate are examined. 
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Zero Spin Rate 


If the spin rate Q is zero, then equation (28) reduces to 


= 0 


(29) 


where 


= (l + j2s^2jSp2 

= [(i ^ ^ ^ + 3 V ^ Sp")]a)^2s, 

Td = j^(j^ - l)^EG(c0r2sj2)2 


2^ 2 
J 


The two roots of equation (29) of greatest magnitude may be calculated by the 
quadratic formula : 


a>Tj 


2 



4TaTd 

V ■ 


(30) 


Since the expression under the radical is never negative, the values of OJ-pj 

are real. For slender rings, the values of and Sp are much less than 

one. Consider the range of mode numbers j such that j^Sp^, and 

j^S ^E_ are much less than one. For this range, equation (30) becomes 
P ^ 



Substitution of equation (31) into equation (26) yields the ratio of the out-of- 
plane deflection w(s,t) to the angle of twist 3(s,t): 


-1 _ Ip F(1 + E(,) 

^ATj AR . 2 

3 + Eg 


(32) 


Equations (31) and (32) describe the out-of-plane inodes which consist primarily 
of twisting motion about the line of centroids accompanied by a small deflection 
perpendicular to the plane of the ring. An example of this type of mode is 
illustrated in figure 1. In this case the deflection of the centroid and the 
rotation about the centroid combine to produce rotation of the cross section 
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about a point which is near the outer edge of the ring. Therefore, the outer 
surface appears to be relatively undistorted, while the inner surface shows a 
large distortion. 

The two roots of equation (29) of least magnitude may be calculated by the 
quadratic formula: 



(33) 


2 2 .22 

For the range of mode numbers such that j Sj- and j Sp Eq are much less 
than one, equation (33) becomes 


(34) 

\lj^ + Eg 


Substitution of equation (34) into equation (27) yields the ratio of the angle 
of twist 6 (s,t) to the out-of-plane deflection w(s,t): 


*^ADj 


R 


1 + Eg \ 

1 + Eg/j2j 


(35) 


Equations (34) and (35) describe the out-of-plane modes which consist primarily 
of deflection perpendicular to the plane of the ring accompanied by a small 
twisting motion about the line of centroids. An example of this type of mode 
is illustrated in figure 2 . 

For mode numbers within the range considered above but significantly 
greater than \J^Qf equations (31) , (32) , (34) , and (35) simplify to the 

following : 


tOTj 




AR 


(1 + Ef:-) 


(36) 


(37) 


to. 


Dj 


= ±j 25 ^ 0)3 


(38) 
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(39) 


'ADj 


= -|(i + Eg) 


where >> E^. 


A 


Mode Number of Zero 

If the mode number j is zero, then equation (28) becomes 


^a^^wo' 


- = 0 


(40) 


where 


T 


a 


2 


Tb = 



+ 



Sr 



2 


The two nonzero roots of equation (40) correspond to a uniform rotation of 
the ring about the line of centroids which has a frequency of oscillation of 


^TO 



The dynamics of this case may be written as 


(41) 


w(s,t) = 0 (42) 

B(s,t) = cos sin co^Qt 

The frequency of this mode is affected by the spin rate. If I^ > then 

the frequency increases with increasing spin rate and has a stiffening effect. 
If I^ < I^, the frequency decreases. 

Figure 3 illustrates this mode in a ring with a rectangular cross section. 
The long axis of the cross section in figure 3 (a) is nominally in the radial 
direction. As the ring spins, the centrifugal force is greater on the outer 
edge of the ring than on the inner edge . This produces a couple which tends to 
force the ring back to its rest condition and thus stiffens the vibration mode. 

The long axis of the cross section in figure 3 (b) is nominally parallel to 
the spin axis. In this case, the couple has an effect opposite to stiffening. 
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The couple may overpower the elastic restoring forces if the frequencies given 
by equation (41) become complex. This occurs when 



and When these conditions are satisfied, the mode is dynamically 

unstable. The lowest spin-rate limit imposed by equation (43) occurs when 
= 0. For this case the limiting spin rate is o)g which is identical to 
that given by equation (14) for static stability. 

For j = 0, the mode frequency given by equation (40) has a double solu- 
tion of zero which implies a solution form different than equations (24) 
and (25) . The dynamics of this case may be written as 

w(s,t) = z + zt (44) 

6(s,t) = 0 

This solution is rigid body translation along the inertial z-direction which 
nominally coincides with the axis of the ring as is shown in figure 4. The 
constant z is the initial displacement and z is the rate of displacement. 


Mode Number of One 

If the mode number j is one, then equation (28) becomes 

Ta“wl^ - = 0 (45) 

where 

Ta = (l + 

= 4(1 + 

= (l + - 2 Sj , 2)(1 + EQ)a3^2sjV + 
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Equation (45) can be factored into 


<“wl - ♦ Vwl" - Vwl - *a) - 0 

where 

= (l + 

Tj, = (l + 2Sr^ + Sz^jd + EG)Wr^Sj2 + (l - Sr^)(Sz^ - 

= (l - 2Sr^ + Sz^)(l + EG)tOr^Sj^^ ~ 

2 2 

If the ring is very slender such that and are very small relative 

to unity, equation (46) may be factored into 






2 


(1 + E^)03^^Sj^ 




0 


(47) 


Two of the roots of equation (47) correspond to the twisting mode which 
has a frequency of oscillation of 


(1 + Eg)Sj2o)^2 ^ (3^2 _ 3 ^ 2^^2 

_____ _ (48) 


2 

Substituting equation (48) into equation (26) with j = 1, assuming that 
and Sp^ are very small relative to unity, and assuming that the spin rate is 
much less than , yield the proportionality constant between the out-of-plane 
deflection w(s,t) and the angle of twist B(s,t): 



C 


-1 

ATI 



(49) 


which is small for slender rings. The spin rate has a small effect on the 
frequency of this mode. 
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For the remaining roots of equation (47) , two different cases may be con- 
sidered. If the spin rate is zero, then the mode frequency has a double solu- 
tion of zero. The dynamics of this case may be written as 

(I) 

(t) 

This solution is a rigid body rotation of the ring about a line lying in the 
XY-plane which nominally coincides with the plane of the ring, as shown in 
figure 5. The constant 0 is the initial rotation angle about the Y-axis 
and 0 is the rate of rotation about the Y-axis. The constant (() is the 
initial rotation angle about the X-axis and (j) is the rate of rotation about 
the X-axis. These angles correspond to the pitch and roll angles of Euler ian 
mechanics so long as their values are small. 

If the spin rate is not zero, then the mode frequency has solutions equal 
to plus and minus the spin rate. The dynamics of this case may be written as: 


3(s,t) = -(0 + 0t) cos 


(t)- 


(cj) + (|)) sin 


w(s,t) = R(0 + 0t) cos [ — j + R((j) + (J)) sin 


w(s, t) 


3(s,t) 



(52) 


(53) 


The first pair of terms of equations (52) and (53) correspond to a constant 
inclination of the ring axis determined by the angles 0 and ([) as shown in 
figure 6. The second pair of terms of equations (52) and (53) describe a gyro- 
scopic precession shown in figure 7. The initial orientation of the ring is 
determined by a rotation about the Y-axis of angle 0p and a rotation about 
the X-axis of angle (j)p. 


Mode Numbers of Two and Greater 

This section examines the vibration characteristics of the out-of-plane 
modes for mode numbers of two and greater when the spin rate is substantially 
less than 0)g and the ring is very slender so that the quantities 
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and j^Sp^ are much less than unity. With these conditions, equa- 
tion (28) can be approximated by 

= 0 (54) 

where 


= 4j^(l + 

= j^(j2 - 1)\(o3^2s^2)2 ^ .2^.2 ^ 


As described in appendix B, the solutions of equation (54) may be 
expressed in the form 


= ^wj ± '^Dj' -"^wj ± “Tj 


(55) 


where 




(56) 


M 2 _ ^2 

‘^Dj - 5 


j + Eg 


(57) 


w, 


^ ^ ^ g ) t ~ 


Tj 


M 2 

r 


(58) 


Since the out-of-plane nodal frequency u^j is much smaller than the vibration 
frequencies, co^j and w.p j , it is ignored in the following discussion. This 
effectively ignores the presence of in equation (54) . 

Equation (58) is the formula for the frequency of vibration of the out-of- 
plane mode which consists primarily of twisting motion about the line of cen- 
troids with small deflection perpendicular to the plane of the ring. The ratio 
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of the out“of-plane deflection to the angle of twist is found by substituting 
equation (58) into equation (26) : 


ATj aR j2 + Eq 


( 59 ) 


The phase velocity relative to the 
ing the modal frequency ((jO,pj/2TT) by the 


rotating ring is obtained by multiply- 
wavelength (27TR/j) : 


Tj 



(60) 


where v^ is the velocity of transverse sound waves. The signs of equa- 
tion (60) are chosen so that the velocities of the traveling waves which move 
around the ring in the same direction as the spin rotation have positive 
values. In this paper, these waves are called progressive. Those waves which 
move in the opposite direction are called regressive. Since the velocities 
given by equation (60) are much larger than the velocity of the ring the 

apparent frequency of oscillation and the phase velocity which would be 
observed from inertial space are effectively the same as those given by equa- 
tions (58) and (60) . 

Equation (57) is the formula for the frequency of vibration of the out-of- 
plane mode which consists primarily of deflection perpendicular to the plane 
of the ring with small rotation about the line of centroids. The ratio of the 
angle of twist to the out-of-plane deflection is found by substituting equa- 
tion (57) into equation (27) : 


1 j2(l + Eq) 


(61) 


The phase velocity relative to the rotating ring is obtained by multiplying 
the modal frequency by the wavelength: 


V 

DD 



(62) 


where Vg is the velocity of longitudinal sound waves . 

The phase velocity relative to inertial space is obtained by adding the 
nominal velocity of a point on the ring to the phase velocity relative to 

the r i ng : 
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V . = + 

SD] 


(j' 


1 )' 


+ E. 




(63) 


The apparent frequency of oscillation observed at a point fixed in inertial 
space is obtained by dividing the apparent phase velocity by the wavelength: 


“SDj ” 




(64) 


For small rotation rates Q, the modal frequencies and phase velocities 
are determined mainly by the ring's material properties (E, G, and p) and 
the geometric parameters (1^,, J, A, and R) . For large rotation rates, the 
modal frequencies and phase velocities approach the following values: 


t^Dj = 

(66) 


(67) 


( 68 ) 


where is significantly larger than 

The restoring forces caused by the tension in the ring are large in com- 
parison to the restoring forces caused by material strain so that given a suf- 
ficiently high rotation rate, the out-of-plane deflection modes of the ring 
approach the solutions for a spinning loop of string (which has = 0) • 

Equation (66) shows that the phase velocities relative to the material approach 
the linear velocity of the edge of the ring in both the progressive and regres- 
sive waves. Equation (67) shows that the phase velocities relative to inertial 
space approach a value twice the ring speed for the progressive wave and a very 
slow velocity for the regressive wave. 

The material of the ring has a velocity of Rfl. The progressive wave 
travels around the ring in advance of the material at a total velocity of 2RS 
as shown in figure 8. The regressive wave shape is almost stationary in space 
with a slight retrograde movement. To an observer, the material of the ring 
would appear to be following a track which is almost fixed in space as shown 
in figure 9. 


Voj = 


SDj 


2RS1 


(j‘ 


1 ) 


j2 + E, 


OJ 




^ iiil 


1) 


SDj 


RH 


+ Eg 
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The following development examines the out-of-plane structural modes in a 
more conventional manner than was done above. Application of trigonometric 
identities to equations (24) and (25) yields the following alternative form for 
the modal equations : 


w(s, t) 


S(s,t) 


oo 

2 (^j ^ 


~ is 

cos + Bj cos sin 


j=2 


Ds 


+ Cj sin ^ sin ^jt + Dj sin 




2 ^^ADj(^j 


is ~ IS. 

cos COS Wgjt + Bj cos ^ sin co^^t 


R 




j=2 

~ is "" is 

+ Cj sin ^ sin cOj^jt + Dj sin cos oOj^jt 


(69) 


(70) 


where the modal frequencies (.Oj-^j are given by equation (57) and the propor- 
tionality constant between the twist angle and the out-of-plane deflection 
C/^Dj is given by equation (61) and where 


A-; = A-; + 


= C- 


- A- 


Bj - Bj Dj 


Dj = Bj + Dj 





(71) 


Therefore, the out-of-plane structural modes may be represented as the sum of 
products of sinusoidal functions of s (mode shapes) and sinusoidal functions 
of time (oscillations) . Since the mode shapes are a function of s only, 
standing waves do not move with respect to the material of the ring; that is, 
the nodes of the standing waves are fixed in the material. 

Consider, for example, any one of the terms of equation (69): 


w(s, t) 


cos 


js 

R 


cos (J^Qjt 


(72) 


Consider an observation point fixed in inertial space and located at a dis- 
tance X along the circumference of the ring measured from the X-axis. The 
deflection of the ring observed at this point is the deflection at 


s = X - R^^t 


(73) 
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Substituting equation (73) into equation (72) yields the motion of an out-of~ 
plane standing wave relative to the inertial coordinate system: 

w(s,t) = Aj cos “(X - RQt) cos t%jt (74) 


In the limit for 0. ^ equation (65) may be substituted into equation (74) . 
Subsequent trigonometric manipulations yield the following representation for 
the modal behavior of the ring at high spin rates as seen by a fixed observer: 


^(X, t) 




n 


cos — X + cos :~(X - 2Rf^t) 


R 




(75) 


This equation consists of two terms. The first is a sinusoidal shape of the 
ring which is fixed in time. The second is a sinusoidal shape which travels 
around the ring at twice the rim speed. Figure 10 illustrates the motion of a 
standing wave described by equation (75) as would be viewed by a fixed observer. 
The ring achieves its maximum deflection in figures 10(a) and 10(g) with inter- 
mediate positions in figures 10(b) through 10(f). The square symbol on these 
figures identifies the reference point fixed in the ring at s = 0. Although 
the crest and troughs (at maximum deflection) alternate when viewed relative to 
a point which is fixed to the ring, they always occur at the same locations in 
inertial space, as can be seen by comparing figures 10(a) and 10(g). 

The analysis of this section has been limited to a range of mode numbers 
for which the quantities r l^^p^ much less than 

unity. Since the effects of transverse shear upon the deformation of the ring 
are neglected in the development of the equations of motion of appendix A, the 
equations derived above may give inaccurate results for mode numbers beyond the 
indicated range. However, for large mode numbers within the range, several of 
the equations can be simplified to the following: 


^Tj “ 

. = 
ATj 



VTj 



0) 


Dj 


±j 




+ 


(76) 


(77) 


(78) 


(79) 
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( 80 ) 


C - " ^ 

^ADj - R 


“SDj = + \/j^Sj.2a)s2 + Q2'j 


VgDj = + \/j^Sj-2vg2 + (RfJ)2 


(81) 

(82) 

(83) 


where >> Eq. 

For high spin rates, equations (82) and (83) become 


0) 


SDj 


2Q 


, -2jQ. 


2„ 2 


SDj 


= “3 


2 Vg 


2R^ ' 


+ 2RT^ 


(84) 


(85) 


IN-PLANE HOMOGENEOUS SOLUTIONS 

The characteristics of the in-plane vibration modes (deflections in the 
radial and tangential directions) are examined in this section. The equations 
of motion of the in-plane modes are given by equations (1) and (2) . In the 
absence of external influences (control forces and gravity) , the distributed 
forces acting in the plane of the ring, W]^ and W 2 , are zero. As in the 
out-of-plane case, the structural perturbation variables are periodic in the 
variable s, that is. 


u (s + 

2TTR) 

= u (s) 

(86) 

V (s + 

2itR) 

= v(s) 

(87) 


for all s. As in all wave equations, there is a family of solutions of equa- 
tions (1) and (2) which satisfy the constraints of equations (86) and (87) . 
Assume that each mode has the solution 


u (s, t) 


cos 


(i" + “ujt) + 


( 88 ) 
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(89) 


v(s,t) = Gj cos + Hj sin 

where the mode number j is a nonnegative integer. 

Substitution of equations (88) and (89) into equations (1) and (2) yields 
the following simultaneous algebraic equations : 


where 


\(j^ - + (j2 _ i)S^2j - (i + 

= [j(ojg2 + 

- (l + 



(90) 

(91) 


The constant Gpj is the ratio between the amplitudes of tangential deflection 
of the centroid and its radial deflection. 


Solving equations (90) and (91) by eliminating Gpj yields a polynomial 
in the in-plane vibration frequency O)^ j : 


*^a^u3 


^ - TkW 

b u j 


- ^c%j‘ 


'^d%j 


Te = 


(92) 


where 


= 1 + (j^ + DS^.^ 

= 4jfiSg2 

= (j^ + 3)^2 + (j2 + l)(wg2 _ + (2 + S2^)(j^ - 

= 4jf2(o)g2 + f^2^ 

Tg = [(j^ - D^S ^OJ + (j2 - 3)W 

^ L z s s J 
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The roots of equation (92) are the vibration frequencies of the in-plane modes. 
The ratios of the amplitudes of the tangential and radial deflections for each 
of these modes are obtained by substituting the values of the frequencies into 
either equation (90) or (91) . 

The remainder of this section first examines the vibration characteristics 
when the spin rate is zero. Then the modes with mode number equal to zero and 
one with nonzero spin rate are examined. Finally, the vibration characteristics 
of the higher numbered modes with nonzero spin rate are examined. 


Zero Spin Rate 

If the spin rate Q is zero, then equation (92) reduces to 

- Tc^uj2 + Tg = 0 (93) 

where 

Tg = 1 + (j2 + l)sj^ 

Tc = [j^ + 1 + (2 + S2^)(j2 - 

Te = j^(j^ - 

The two roots of equation (93) of greatest magnitude may be calculated by the 
quadratic formula: 


^Cj 


2 




(94) 


Since the expression under the radical is never negative, the values of co^j 
are real. For slender rings, S^ is much less than one. Consider the range 

o p 

of mode numbers j such that j is much less than one. For this range, 

equation (94) becomes 


j = ± n/: 


+ 1 CO, 


(95) 


Substitution of equation (95) into equation (91) yields the ratio of the 
tangential deflection v(s,t) to the radial deflection u(s,t): 


^FCj ^ 


(96) 
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Equations (95) and (96) describe the in-plane modes which consist primarily of 
compression and stretching of the line of centroids with little bending. An 
example of this type of mode is illustrated in figure 11. 

The two roots of equation (93) of least magnitude may be calculated by the 
quadratic formula 


“Bj 


2 



(97) 


For the range of mode numbers such that is much less than one/ equa- 

tion (97) becomes 


CO 


Bj 


= + 


j(j^ - 1) 
+ 1 


^z^s 


(98) 


Substitution of equation (98) into equation (91) yields the ratio of the radial 
deflection u(s,t) to the tangential deflection v(s,t): 



(99) 


Equations (98) and (99) describe the in-plane modes which consist primarily of 
bending of the line of centroids with little compression and stretching. An 
example of this type of mode is illustrated in figure 12. 

For mode mombers within the range considered above but significantly 
greater than one, equations (95) and (98) reduce to the following: 



(^cj = ±:i'^s 

(100) 


<^Bj = 

(101) 

where 

» 1. 



If the mode number j 


^a“uO^ - ^c%0^ 


Mode Number of Zero 
is zero, then equation (92) becomes 

0 
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(102) 


,1 


I 


where 


Ta = (l + s/) 

= (l + ^ ^3 _ 

p 

For slender rings, the constant may be ignored. The two nonzero roots of 

equation (102) correspond to a radial oscillation of the ring with no bending 
which has a frequency of oscillation of 

“cO "" - \j“s^ 


Substitution of equation (103) into equation (91) yields the ratio of tangential 
deflection to radial deflection: 


%co 


2^^ 


/ojg2 + 31^2 


(104) 


The relative magnitudes and phasing of the radial and tangential terms are such 
that the angular momentum of the ring is constant. 

The radial oscillation mode is illustrated in figure 13. ^ For this figure 
the parameter F^q (eq. (88)) is zero. Figure 13(a) shows the ring at its 
initial condition (t = 0) . The radius is increased by the amount Eqq. The 
tangential velocity of the ring is the sum of the nominal velocity, R^, and 
the time derivative of the variable v(s,t). Since the radius of the ring is 
larger than the nominal radius at t = 0, the ring rotates at a speed less than 
nominal. Figure 13(b) shows the ring after one-quarter cycle of oscillation 
(t = tt/2Dqq) . The radius of the ring equals the nominal radius and is decreas- 
ing. Because the ring was rotating at a slower than nominal rate previously, 
it has retarded by an amount ^fcO^CO* This retardation is shown by the loca- 
tion of the mass element at s = 0 (denoted by the dark square fixed to the 
ring) relative to its nominal location (denoted by the light square) . The ring 
is now rotating at its nominal rate and is decelerating. Figure 13(c) shows 
the ring after one-half cycle of oscillation . (t = e/o)^q) . The radius of the 
ring is now less than nominal by the amount E(-;q; the ring is rotating at a rate 
which is greater than nominal and the previous retardation has been canceled. 
Figure 13 (d) shows the ring after three-quarters cycle of oscillation 
(t = 37T/200qq) . The radius of the ring is now equal to its nominal value? the 
ring is rotating at its nominal rate and is decelerating. The indicated mass 
element has progressed ahead of its nominal circumferential location. 

The double root at zero in equation (102) implies a solution which is not 
of the form of equations (88) and (89) . For this case, the following form may 
be used: 
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u (s, t) = R 


(105) 


v(s,t) = t/; + 


(106) 


Substituting equations (105) and (106) into equations (1) and (2) yields 


(l + 


(107) 


ip = Arbitrary value 


(108) 


The parameter specifies an arbitrary angular displacement of the ring about 

its spin axis. The parameter specifies a constant increment in the spin 

rate of the ring; that is, the actual spin rate is equal to Q + The 

parameter R is the change in the radius of the ring caused by the change in 
centripetal force acting on the ring. 


Mode Number of One 

If the mode number j is one, then equation (92) may be written in the 
factored form 


(Wui (l + 


+ 203 ^^ ) = 0 


(109) 


The roots of the quadratic term of equation (109) are 


± + (l + + 20 ) 3 ^) 


1 + 2S. 


( 110 ) 


Substitution of equation (110) in equation (91) yields the ratio of the tangen- 
tial deflection to the radial deflection: 


^FCl “ ^ 


( 111 ) 


Assuming that the spin rate is substantially less than (jO^ and that the param- 
eter is much less than unity reduces equation (110) to 


= is/Ff 


( 112 ) 
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This solution is a compression mode in which the ring maintains its circu- 
lar shape but moves in such a way as to keep its center of mass fixed in space. 
This mode for the negative value of equation (112) and Fq-|^ (eq. (88) ) equal 
to zero is illustrated in figure 14. At time zero (fig. 14(a)), the shape of 
the ring is shifted in the positive x-direction while most of the mass is dis- 
placed in the negative x-direction. After one-half cycle of oscillation 
(fig. 14(c)), the ring's shape has rotated counterclockwise to the negative 
x-direction and the mass has accumulated in the positive x-direction. For the 
positive value of equation (112), the motion is similar except that the ring's 
shape rotates clockwise. The sum of these two modes produces the standing-wave 
motion of figure 15 where the ring is oscillating from side to side. Fig- 
ures 14 and 15 which illustrate this mode for zero spin rate are representative 
of the motion of this mode for nonzero spin rate because the frequency of oscil- 
lation of the mode is so much greater than 

The double root in equation (109) implies a solution which is not of the 
form of equations (88) and (89) . For this case, the following form may be 
used : 


u (s, t) 
v(s, t) 


(x + xt) 


(y + yt) 


cos 


cos 


(i * * 

(I " - 


(y + yt) 
(x + xt) 


sin 


sin 


(I * 

(t ^ 


(113) 

(114) 


These solutions represent rigid body translations and rates in the plane of the 
ring. They are illustrated in figure 16. 


Mode Numbers of Two and Greater 

This section examines the vibration characteristics of the in-plane modes 
for mode numbers of two and greater when the spin rate is substantially less 
than (jjg and the ring is very slender so that is much less than unity. 

In addition, assume that the roots of equation (92) are each much greater than 
4jQS2^. With these assumptions, equation (92) can be approximated by 


where 

Tc = (j2 + 

2 

T^ = 4jf20)g 

+ (j^ - 


(115) 
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As described in appendix B, the solutions of equation (115) may be expressed in 
the form 


CO. 




= Uj 


± (jO 


'Bj ' 


-Uj 


± CO, 




(116) 


where 


U. = — ^ (117) 

^ + 1 

“Bj - 

“Cj ' \/jB + 1 Wg (119) 


Define the values of which are solutions of equation (115) as 

follows : 


‘^BRj 

‘^BPj 

“CRj 

‘^CPj 


<^Bj + 


■^B3 + "j 


‘"cj - 


-OJcj - Uj 


(120) 

( 121 ) 

( 122 ) 

(123) 


Equations (120) and (121) with equations (117) and (118) are the formulas 
for the frequencies of vibration of the in— plane modes which consist primarily 
of bending of the line of centroids with little compression and stretching. 

The ratios of the amplitudes of the tangential movement v(s,t) to the radial 
deflection u(s,t) are found by substituting equations (120) and (121) into 
equation (91) : 


^FBRj 



(124) 


%BPj 



( 125 ) 


where Gpgj^j and Gpgpj are the ratios corresponding to the frequencies 
and Wgp j , respectively. 
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The phase velocities relative to the rotating ring are obtained by multi- 
plying the modal frequencies (00gj^j/2TT and 00gpj/27T) by the wavelength (2'TTR/j): 


V- 


(j2 - 1) \j + DS^^Vg^ + (Rfl)2 + 


BRj 


BPj 


3 ^ -HI 


(j^ - 1) \J(j^ -*• + (R^^)2 - 2Rf2 


j2 + 1 


(126) 


(127) 


where v^ is the velocity of longitudinal sound waves. The signs of equa- 
tions (126) and (127) are chosen so that the velocities of the traveling waves 
which move around the ring in the same direction as the spin rotation (progres- 
sive waves) have positive values. 

Consider the solution of equations (88) and (89) where Fj = 0, 03^ j = 

and ^ << S 2 ^s. Substitution of equation (121) into equations (88) and (89) 
and using equations (117), (118), and (125) yield the following: 


u(s,t) = Ej cos 


V (s 


\f- (‘"Bj - 

,t) = sin - (OJgj - Uj)t] 


where 


""j = ^ 

+ 1 


,, . j (j^ - 1 ) ^ 

COBj — 




j2 + 1 


(128) 


(129) 


(130) 


(131) 


Equations (128) and (129) describe a sinusoidal traveling wave which progresses 
ahead of the nominal material motion with the relative velocity 


"^BPj ^Bj 

where 


VBi = S V 

bd z s 


\|j^ -H 1 


V • = 


' + 1 




(132) 


(133) 


(134) 
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This mode consists mainly of bending of the line of centroids with little com- 
pression and stretching. Since the spin rate is sxibstantially less than 
the vibration frequencies and phase velocities are primarily functions of the 
geometric and material properties of the ring (E, p, A, 1 ^, and R) and the 
mode numbers. Figure 17 illustrates this mode for j = 2. 

Similarly, consider the solution of equations (88) and (89) where Fj = 0, 
0)uj = j ' ^ ^z^s- Substitution of equation (120) into equations (88) 

and (89) and using equations (117), (118), and (124) yield the following: 


u (s, t) 


v(s,t) 



(0)Bj 




+ Uj)^ 
+ Uj)^ 


(135) 


(136) 


where Uj and OJgj are given by equations (130) and (131) . 

Equations (135) and (136) describe a sinusoidal traveling wave which 
regresses behind the nominal material motion with the relative velocity 


V- 


BRj 


Bj 


+ V ■ 


(137) 


where v^j and Vj are given by equations (133) and (134) . Figure 18 shows 
this mode for j = 2. 

Since C 03 j is substantially larger than the spin rate for mode numbers of 

two and greater, the apparent phase velocities, v • and v . , and the 

SBPj SBR]i 

apparent vibration frequencies, and , which would be seen by a 

stationary observer are essentially equal to j , Vgg j , Wgp j , and j , 

respectively. 

Given that the ring has both progressive and regressive in-plane bending 
modes of equal amplitude, the resultant motion would have a standing-wave 
character. Adding the two traveling waves described by equations (128) , (129) , 

(135) , and (136) and applying trigonometric identities yield the following 
equations for the in-plane bending-mode standing waves: 


u(s,t) = Ej cos + Ujtj cos Wgjt (138) 

v(s,t) = sin t Aijtj cos ^Bjt (139) 

where Uj and U)gj are given by equations (130) and (131), Ej = 

^ « S^cOg. Following the above analysis, using Ej = 0 yields equations 
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similar to equations (138) and (139) except for interchanges in the sine and 
cosine operators. Therefore, the in-plane structural bending modes may be 
represented by the sum of products of sinusoidal functions of time (oscilla- 
tions) and sinusoidal functions of distance along the ring s and time (mode 
shapes). Since the mode shapes are functions of time as well as s, the stand- 
ing waves move with respect to the material of the ring. The nodes of the 
standing waves may be Said to move along the ring at some nodal velocity. The 
nodal velocity of the bending standing waves is equal to Vj (eq. (134)) in 
the direction opposite to the rotation of the ring. In other words, the modes 
regress at the velocity Vj relative to the material. At time equal to zero, 
the deformation variables as seen by a fixed observer can be plotted as shown in 
the upper two curves of figure 19. The radial displacement u is greater than 
the tangential displacement v by the factor j . After one-quarter cycle of 
oscillation, the plots for both variables collapse to zero, corresponding to an 
undeformed ring (not shown). After one-half cycle of oscillation (t = , 

the variables u and v attain the plots shown in the lower two curves of 
figure 19. The ring and the s coordinate system have moved by the amount 
I^t. The nodes of the mode shape have lagged behind by an amount Vjt- Fig- 
ure 20 shows the ring with a bending-mode standing wave for j = 2. After one 
full cycle of oscillation (t = 27r/(j0gj) , the ring has rotated through an angle 
of f2t and the structural mode shape has lagged behind by an angle of (Vj/R)t. 


For higher spin rates, the character of the bending modes changes. The 
above analysis considered low spin rates for which The following 

<< Q « 


analysis examines higher spin rates for which 


^z^s 


s • 


For this range 


of spin rates, the above development applies except for the formulas for o) 


and V 


Bj 


which become 


'Bj 


(jO 


Bj 



- 1 ) 
+ 1 


a 


V 


Bj 



(140) 


(141) 


where << 0 . « co^. The frequencies and phase velocities of the bending 

modes are now functions of the tension in the ring caused by centrifugal force. 

The progressive bending -mode traveling wave described by equations (128) 
and (129) for high spin rate and j = 2 is illustrated in figure 21. This 
figure is similar to the figure for low spin rate (fig. 17) except for the 
speed of rotation of the ring's shape and the movement of the ring's material 
relative to the ring's shape. For low spin rates, the ring's shape rotated at 
a relatively slow rate which was significantly faster than the nominal rotation 
rate of the ring. For high spin rates, the ring's shape rotates at a relatively 
high rate which is at most twice as fast as the nominal rotation rate of the 
ring. For j = 2 (shown in fig. 21) , the relative phase velocity Vgpj is a 
fraction of the material velocity. The ring may thus appear to have a fixed 
distortion and to be rotating as a rigid body. This illusion diminishes with 
increasing mode number. 
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The regressive bending-mode traveling wave described by equations (135) 
and (136) for high spin rate and j = 2 is illustrated in figure 22. This 
figure is similar to figure 18 for low spin rate except that now the speed of 
the wave shape relative to the material is equal to the speed of the material 
in the opposite direction. Therefore, the shape of the ring is fixed in space 
and the material moves around this fixed contour . 

The phase velocities of these modes relative to inertial space, 
and Vg 3 j^j , are obtained by adding the nominal velocity of the ring to 

the phase velocities relative to the ring, Vgpj and Vgp j , given by equa- 
tions (126) and (127) : 


^SBRj 


= -(j 


^2 


1 ) 


S ^v ^ 
2R^^ 


(142) 


^SBPj ~ 



(143) 


The apparent frequencies of oscillation, ^sBPj ^SBR j ' observed from 

inertial space are obtained by dividing the apparent phase velocities by the 
wavelength: 


‘^SBRj 


j (j 


.2 



(144) 


“^SBPj 



(145) 


Given that the ring has both progressive and regressive in-plane bending 
modes of equal amplitude and << << oOg, equations (138) and (139) apply 

except that Wgj is given by equation (140) - Figure 23 shows the deformation 
variable as seen by a fixed observer at time equal to zero and after one-half 
cycle of oscillation. Using equation (140) , the time for one-half cycle of 
oscillation is 




+ 1 
- 1 ) 


TT 


(146) 


The material of the ring has moved a distance equal to 


‘ 4 - 1 

Rflt = R7T 

j(r - 1 ) 


(147) 
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Using equation (134) , the distance that the mode of the standing wave has 
regressed relative to the material is 


V.t = RTT (148) 

j(j - 1) 

Subtracting equation (148) from (147) gives the spatial displacement of the 
mode of the standing wave: 

Rf2t - Vjt = J fT (149) 

which is equal to one-half of the wavelength of the structural deformation. 
Therefore, the peaks and valleys of the structural vibrations occur at the same 
locations in space and do not interchange as was the case for very low spin 
rates shown in figure 19. Figure 24 shows the ring with a bending-mode stand- 
ing wave for j - 2, 

Equations (122) and (123) with equations (117) and (119) are the formulas 
for the frequencies of vibration of the in-plane modes which consist primarily 
of compression and stretching of the line of centroids with little bending. 

The ratios of the amplitudes of the tangential movement v(s,t) to the radial 
deflection u(s,t) are found by substituting equations (122) and (123) into 
equation (91) : 


ScRj = (150) 

G^CPj = -3 (151) 

where *^FCPj ratios corresponding to the frequencies 

and j < respectively. 

The phase velocities relative to the rotating ring are obtained by multi- 
plying the modal frequencies (03(^j^j/2TT and oo^pj/2TT) by the wavelength (2TTR/j): 


1 2R^ 

•CRj ~\l j2 .2^^ 


, 1 2Rf^ 

1 \ + 1 


(152) 

(153) 


Since these velocities are much larger than the velocity of the ring Rf^, the 
apparent frequencies of oscillation, ^SCRj ' apparent phase 
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velocities, v„-d • , that would be observed from inertial space are 

effectively the same as those given by equations (122) , (123) , (152) , and (153) . 

Consider the solution of equations (88) and (89) where Fj = 0, 

OJ^j = c^Qpj, and Q << 0)^. Substitution of equation (123) into equations (88) 
and (89) and using equations (117), (119), and (151) yield the following: 


u(s,t) = Ej cos 


f- 


v(s,t) = “jEj sin 


[x'“ (‘"Cj + 


where 


u . = 


2j 


' - 1 




“cj ■ 


y + 1 Ui 


(154) 


(155) 


(156) 


(157) 


Equations (154) and (155) describe a sinusoidal traveling wave which progresses 
ahead of the nominal material motion with the relative velocity 


^CPj = ^Cj + 


where 


^cj = f ^ j2 


v ■ = 






(158) 


(159) 


(160) 


This mode consists mainly of tension and compression distortions with very 
little bending of the ring’s material. The phase velocity is primarily a func- 
tion of the material and geometric properties of the ring. 

Similarly, consider the solution of equations (88) and (89) where Fj = 0, 
j / and << oOg. Substitution of equation (122) into equations (88) 

and (89) and using equations (117), (119), and (150) yield the following: 

u(s,t) = Ej cos + (0)Qj - (161) 
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(162) 


v(s,t) 


-jEj sin 



<^Cj 



where Uj and oo^j are given by equations (156) and (157) . Equations (161) 

and (162) describe a sinusoidal traveling wave which regresses behind the 
nominal material motion with the relative velocity 


"CRj 


= V, 


Cj 


(163) 


where v^j and Vj are given by equations (159) and (160) . 

Given that the ring has both progressive and regressive in-plane compres- 
sion modes of equal amplitude, the resultant motion would have a standing-wave 
character. Adding the two traveling waves described by equations (154) , (155) , 

(161) , and (162) and applying trigonometric identities yield the following equa- 
tions for the compression-mode standing waves: 


u (s, t) 
v(s, t) 


E-: cos 

P - u.t) 

cos 

“cj 

3 

V R 3 / 


— T- Sin 

- u.t) 

cos 

«Cj 

3 

\ R 3 / 


(164) 

(165) 


where Uj and U)^j are given by equations (156) and (157), Ej = 2E j , and 

Q « 0 ) 3 . Since the mode shapes are functions of time as well as s, the stand- 
ing waves move with respect to the material of the ring. The nodal velocity of 
the compression standing waves is equal to Vj (eq. (160) ) in the direction of 
rotation of the ring. In other words, the modes progress at the velocity Vj 
relative to the material. 


At time equal to zero, the deformation variables as seen by a fixed 
observer can be plotted as shown in the upper two curves of figure 25. The 
tangential displacement v is greater than the radial displacement u by the 
factor j. After one-quarter cycle of oscillation, the plots for both variables 
collapse to zero, corresponding to an undeformed ring (not shown) . After one- 
half cycle of oscillation (t = Tt/oj^j), the variables u and v attain the 
plots shown in the lower two curves of figure 25. The ring and the s coordi- 
nate system have moved by the amount R^t. The mode shape has moved an addi- 
tional amount ^jt. Figure 26 shows the ring which has a compression-mode 
standing wave for j = 3. After one full cycle of oscillation (t - 2Tr/a)^j), 

the ring has rotated through an angle of f^t and the structural mode shape has 

moved through an additional angle of (Vj/R)t. 

The analysis of this section has been limited to a range of mode numbers 

for which is much less than unity. Since the effects of transverse 

shear upon the deformation of the ring are neglected in the development of the 
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equations of motion of appendix A, the equations derived above may give inaccu- 
rate results for mode numbers beyond the indicated range. However, for large 
mode numbers within the range, several of the equations can be simplified to 
the following : 

u . = ~ 0. 




(166) 

% = j + ^2 




(167) 

Wcj = j^s 




(168) 

^BRj = ■ 




(169) 

'^BPj = 




(170) 

^CRj ^ 




(171) 

^CPj "" ^s 




(172) 

2 

where j >> 1. For low spin rates, the formulas 

for the 

bending modes 

become 

"^Bj = 




(173) 

. = -is v 
BRj z s 




(174) 

^BPj ^ j^z^s 




(175) 

where >> 1 and Q « 5^03^. For higher spin rates, the formulas for 

bending modes become 

the 

t^Bj = 




(176) 

^BRj = 




(177) 

VBpj = 




(178) 
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(179) 


where 


SBRj 

2Rfi 

^SBPj 

= 2Rf2 


-;2g 2 (jq 2 

‘^SBRj 

20. 

“SBPj 

= -2t0 

» 1 

and << 


(180) 

(181) 

(182) 


EXAMPLE 


Solutions for an example ring were computed using the formulas of the 
previous sections. The example ring has a radius of 360 m and a 1 . 4- cm- square 
cross section. The mass density of the material is 1.8 g/cm^. Young's modulus 
is 280 GPa, and the shearing modulus is 110 GPa. The other physical character- 
istics of the ring are given in table I. 


The highest spin rate considered for this example is 1 rad/sec. This 
spin rate produces a steady-state strain of the line of centroids of 

= 8.338 X 10“"^ as computed by equation (13). The maximum steady-state 
stress is Cfrnax “ 0.2335 GPa as computed by equation (15). The validity of 
the approximations of equations (A50) through (A53) used in deriving the equa- 
tions of motion in appendix A have been checked. 


The out-of-plane frequencies, ^ j / j / and Wp j ^ were computed from 


equation (28) using the procedure of appendix B. The nodal frequencies u, 


are plotted in figure 27. The vibration frequencies of the torsional modes 


wj 


03 | 


iipj and of the out-of-plane deflection modes co^j are plotted in figure 28. 
The ratios of the out-of-plane deflection to the twist angle for the torsional 


modes 


.“1 

^ATj 


were computed by equation (26) . The ratios of the twist angle to 


the out-of -plane deflection for the deflection modes were computed by 

equation (27). These ratios for = 1 rad/sec are plotted in figure 29. The 
phase velocities for the out-of-plane modes are plotted in figures 30 and 31. 


The in-plane frequencies, u j , j ' were computed from equa- 

tion (115) using the procedure of appendix B. The nodal frequencies Uj are 
plotted in figure 32 and the vibration frequencies, oo^j and j ' plotted 

in figure 33. The ratios of the amplitudes of the tangential motions to the 
amplitudes of the radial deflections were computed from equation (90) for com- 
pression modes and from equation (91) for bending modes. These ratios for 

= 1 rad/sec are plotted in figure 34. The nodal velocities for the in-plane 
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modes are plotted in figure 35. The phase velocities for the in-plane vibra- 
tion modes are plotted in figures 36 and 37. These velocities are relative to 
the nodal velocities. 

Note that the quantities plotted in figures 27 through 37 are functions of 
only integer values of the mode number j. The curves have been faired to 
improve their readability. 

For the example problem being considered here, the out-of-plane nodal 
frequencies u^j are small and do not significantly contribute to the solution 
of the problem. The in-plane nodal frequencies substantially contribute to the 
vibrational behavior of the ring for low-numbered modes. 

The frequencies of the torsional modes (jOipj and the compression modes 0)^^ 
are roughly proportional to mode number. Since the cross section of the ring is 
square (I^ = I^) , the frequencies of the torsional modes are independent of spin 
rate. The phase velocities of the torsional modes Vrpj are approximately equal 
to the speed of transverse sound waves in the material of the ring. The phase 
velocities of the compression modes v^j are approximately equal to the speed 
of longitudinal sound waves. 

The frequencies of the out-of-plane deflection modes and of the 

in-plane bending modes C 03 j are functions of both mode number and spin rate. 
Given low spin rates or high mode numbers, the restoring moments caused by bend- 
ing are larger than the restoring moments caused by tension in the ring. Under 
these conditions, the ring vibrates much like a linear beam. The vibration 
frequencies are proportional to mode number squared and are independent of 

spin rate. The phase velocities are proportional to mode number. For high spin 
rates and low mode niombers , the restoring moments caused by tension dominate so 
that the ring vibrates much like a flexible string held under tension. The 
vibration frequencies are proportional to mode number and to spin rate. The 
phase velocities approach the velocity of the ring rQ. 

The out-of-plane deflection w(s,t) and the twist angle 3(s,t) vibrate 
in phase. For the deflection mode, the ratio of the twist angle to the out-of- 
plane deflection ranges from -1/R to - (1 + ‘Sq) /R rad/m. For the torsional 
mode, there is very little out-of-plane deflection. 

The tangential movement v(s,t) and the radial deflection u(s,t) vibrate 
in quadrature. For the in-plane bending mode, the ratio between the amplitudes 
of the tangential movement and the radial deflection is approximately equal to 
the inverse of the mode number. For the compression mode, this ratio is 
approximately equal to the mode number. 


CONCLUSIONS 

Small -amplitude dynamics of a steadily rotating slender ring have been 
examined. Derived were linear partial differential equations which describe 
flexural and extensional motions in the plane of the ring (in-plane bending and 


42 



compression modes ) , flexural motions perpendicular to the plane of the ring 
(out-of-plane bending modes ) , and twisting motions about the centroid of the 
ring cross section (torsional modes) . The following conclusions were drawn 
from the investigation: 

1. A maximum spin rate, beyond which the ring is statically unstable, has 
been found. Given an increase in the radius of the ring, if the tension change 
required to balance the centrifugal force is greater than the increase in ten- 
sion caused by stretching, then the radius will grow without bound. This 
instability occurs when the spin rate approaches the speed of longitudinal sound 
waves divided by the radius of the unstressed ring. For practical materials 
such as steel, the ring will fail in tension well before this spin rate can be 
achieved. 

2. The spin rate can have an effect upon the vibration frequency of the 
torsional modes. If the cross section is such that the axis of greatest moment 
of inertia is parallel to the radius vector from the center of the ring, then 
the torsional modes are stiffened by increases in spin rate. But, if the axis 
of greatest moment of inertia is perpendicular to the radius vector, then the 
torsional modes are softened by increases in spin rate. The torsional mode with 
mode numbers of zero becomes dynamically unstable at a spin rate which is near 
the limit for static stability. 

3. The spin rate affects the nature of the bending modes of the ring. 

Given low spin rates or high mode numbers , the ring behaves like a vibrating 
beam. For high spin rates and low mode numbers, the ring vibrates like a 
flexible string under tension. 

4. Standing waves (mode shapes) move with respect to the ring. For the 
out-of-plane modes, this nodal motion is insignificant. This motion is signifi- 
cant for in-plane bending modes with low mode numbers. 

5. Some modes of a spinning ring would appear to a fixed observer to be 
surface movement along a fixed, warped, contour. Such modes might pose a prob- 
lem in the design of an active control system because of difficulties in obtain- 
ing the rates of change of the perturbation variables. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 5, 1980 
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APPENDIX A 


DERIVATION OF EQUATIONS OF MOTION 

The linearized equations of motion for a steadily rotating slender ring 
are derived in this appendix. The deformations of the ring are defined by the 
three-dimensional displacement of the centroid of the ring's cross section and 
the angle of twist of the ring about the line of centroids- The stresses 
within the ring and the linear and angular accelerations of each mass element 
of the ring are derived as functions of the ring deformations. Deformations 
caused by shear are not included in this development. Euler's equation and 
Newton's law are used to obtain a linearized system of partial differential 
equations which define the dynamics of the perturbation variables with respect 
to time and location around the circumference of the ring. 


Geometry of Ring 

The geometry of the ring is shown in figure 38. The ring is a slender 
piece of material formed into the shape of a hoop with radius R. The ring 
lies in the plane of the X and Y inertial axes. Choose the point P, which 
happens to be at the intersection of the ring with the X-axis, as a reference 
point; then any point Q on the ring may be identified by the variable s which 
is set equal to the length of the arc PQ. Given a steady spin rate the 

wing expands to a new radius, R + AR. The angle between line OP and the X-axis 
is equal to J^t- Since the ring, being uniform, expands equally in all direc- 
tions, the angle POQ does not change and the length of the arc PQ becomes 
(1 + AR/R)s. 

The displacement of the centroid of the ring caused by structural deforma- 
tions is shown in figure 39. Given the perturbation variables u, v, and w, 
the point Q moves a distance v along the arc of the undeformed ring to Q^^, 
a distance u along a radius to Q 2 r and a distance w parallel to the Z-axis 

to Q 3 . 


The vector which locates the point Q^ in inertial space may be written 
as 


rc(s,t) 


(R + Ar + u) ^cos (I + i 

* (I * 


(Al) 


where u, v, and w are functions of s and t and i, j, and k are unit 
direction vectors parallel to the X-, Y-, and Z-axes, respectively. Given that 
u, v, and w are very small, equation (Al) may be written in terms of cylin- 
drical coordinates : 
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r^(s,t) = (R + AR + u)e^ + VS0+ we^. (A2) 

where the cylindrical coordinate system is defined at the point Q as shown in 
figure 39 . The velocity and acceleration are written as follows: 

rj,(s,t) = (u - f 2 v)ej. + + ft(R + Ar + u)J e0 + we^ (A 3 ) 

^c(s,t) = & - 2 f 2 v - (R + Ar + u)^ 6j. + [^ + 2 f^u - + we2 

(A 4 ) 

where a dot over a variable denotes differentiation with respect to time. 


Forces and Moments 

A curvilinear coordinate system is developed to aid in deriving the stress 
and strain of the deformed ring. Let the vector 62 be defined tangent to the 
deformed line of centroids of the ring as shown in figure 40 . Then, 


®2 = 



ee 


+ 


R 

R + Ar 




(A 5 ) 


where the prime denotes differentiation with respect to s and where |r^* | 
has been approximated by 


R + Ar 
R 


+ v' 


+ 


u 

R 


Let the vector e^^ be^defined perpendicular to the vector e2 and the vector 
obtained by rotating e^ about e2 through Sr the angle of twist of the 
cross section (see fig. 40 ) - Then 


^1 


/S. 





e^ - 


R 


R + Ar 


(u- - 




(A6) 


/S /V 

Let the vector 03 be defined perpendicular to the vectors and e2 so 

that the triad (61,62^33) forms an orthogonal coordinate system: 





X e. 


+ 3e^ - 


r + Ar 


(A 7 ) 


45 


I 



APPENDIX A 


Equations (A 5 ) through (A 7 ) define vectors which have lengths of unity 
plus terms which are products of the perturbation variables; that is, e^^, 
and are unit vectors to the first order. Since the perturbation variables 

are assumed to be small quantities and nonlinear terms are ignored throughout 
this development of linearized equations, the vectors e^, 62/ and 03 are 

treated as unit direction vectors . 

Assume that the cross section of the ring is perpendicular to 63 and the 
principal axes coincide with ej and 03 . Let the components of the forces 
acting on the cross section of the ring be aligned to the (e]^,e2/e3) 

triad- The component F2 is the tensile force in the ring and Fj and F3 
are shear forces. Let any externally applied influences be represented by the 
distributed forces (force per distance) W]^, W2, and W3 aligned to the 

(61,02/03) triad. Given a short segment of the ring of length ds, the sum of 
forces acting on the segment is 


dF = 


d ( ^ 
3s 




+ F 363 


) + 


^2^2 


^ 3^3 


ds 


(A8) 


Ignoring products of u, v, w, and 3 with F]_ * , F2*/ and F3 * and using 

the relations of equations (A 5 ) through (A 7 ) yield an expression for the spatial 
derivative of the total force acting on the ring segment in terms of cylindrical 
coordinates : 


+ Fn ' “ ^ F2 + ■ " 

ds 1 1 1 R 2 R + AR 


+ <W2 + F2- + I Fi - 


[('^' - 1)(«2 + i ^ 1 ) + - T-)^2] + ew3 + 6-F3ja^ 

(u- - - i ^2) + + WW3 + w"F, 


+ i 


[”3 + F3' + + «"^2) - - B-F^Ji 


(A 9 ) 


The strain in a fiber of the ring located at + ^f^3 relative to 

the centroid of the cross section is given by 


r^' (X r, Z f ,u,v,w, 3 ,Ar) 

e(Xf,Zf) = .ri 1 

V* * 


r^’ (x^ , , 0 , 0 , 0 , 0 , 0 ) 


(AlO) 


where 


= r 


+ + 
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An expression for the distributed strain perpendicular to the ring's cross 
section is obtained by substituting equations (A2 ) , (A6 ) , and (A7) into equa- 

tion (AlO) and using small-value approximations for the perturbation variables: 



(All) 


The stress is assiomed to obey Hooke's law: 


a(x^,z^) = E e(x^,z^) 


(A12) 


where E is the modulus of elasticity of the material. It is assumed that 
there is no intrinsic damping in the material of the ring; that is, deformations 
are conservative. Such damping would appear in equation (A12) as an additional 
term which could be a function of e. The total tensile force F 2 is obtained 
by integrating equation (A12) over the cross section: 


F 


2 



e dA 


(A13) 


An infinite series expansion of the denominator of equation (All) is used, all 
moments of the cross section above the second order are discarded, and cross 
products of the cross section are ignored to give an expression for F 2 : 



(A14) 


The second term of equation (A14) arises from the fact that the ring is curved 
and the elastic center does not coincide with the centroid. For very slender 
rings, this term may be ignored. 

The moments caused by the tensile forces acting on the cross section are 
obtained by integrating the cross product of the vector from the centroid to 
each fiber with the stress vector of the fiber: 

(^f®l ^f^3^ ^ E€e2 dA (A15) 

Performing this integration, ignoring the moments of the cross section above 
the second order, and ignoring cross products of the cross section yield 
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where 


M = 


14 = El f w" - — ^ 

1 ^^r\R + AR r; 


M3 = -EIz 


PVfu" - 4. H 4. v.j] 

[r + ar\ r / rVr r Jj 


(A16) 


(A17) 


(A18) 


/s. 

where are the second moments of the cross section about 03 

and ej_, respectively. The second term of equation (A18) arises from the fact 
that the ring is curved and the elastic center does not coincide with the 
centroid. 


The angle of shear strain caused by twisting of the ring y can be deter- 
mined for any fiber in the ring by computing the angle between the fiber and the 
line of centroids (ref. 15) : 


rf * X r^' (0,0) | 

rf ■ (Xf,Zf) I • |r^* (0,0) 


(A19) 


where ^f^l ^f^3* Noting that r^* (0,0) = r^* and using equa- 

tion (A5) yield a simpler form for equation (A19) : 


y(x^, z^) 




+ 

|r^' (x^,z^) 


(A20) 


Using equations (A2) and (A5) through (A7) and assuming that there is no warping 
of the cross section yield an expression for shear strain acting in the ring's 
cross section: 


^ (A 21 ) 

Because of warping, this equation is exact only for circular cross sections. 

The shearing stress is obtained by multiplying equation (A21) by the shear- 
ing modulus of elasticity G: 


t(x^,z^) g y(x^,z^) 


(A22) 
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As above, it is assumed that there is no intrinsic damping in the material of 
the ring. Such damping could appear in equation (A22) as an additional term 
which would be a function of y. 

The line of action of the shearing stress is perpendicular to the radius 
vector from the centroid to the (x£,Zf) location. The moment produced by twist- 
ing is the product of shearing stress and the radius, integrated over the cross 
section: 


M 


2 






T (Xf ,Zf ) 


dA 


(A23) 


Using equations (A21) and (A22) in equation (A23) yields an expression for the 
torsional moment in the ring expressed in terms of the perturbation variables: 


M 


2 


GJ 


R 

R + AR 


+ 


w* \ 
R + AR/ 


(A24) 


where J is the torsional constant for the cross section. For rings which have 
circular cross sections, J becomes the polar moment of the cross section. For 
rings which have noncircular cross sections, J is the torsional constant which 
accounts for the warping of the cross section. A table of torsional constants 
J for typical cross sections is given in reference 10. 

The relationships between the internal forces and moments of the ring are 
now derived according to the procedure of reference 9. Consider a short seg- 
ment of the ring having length ds . Application of Euler's equation to this 
segment yields the relationship between the angular rate of change of the seg- 
ment and the forces and moments acting upon this segment, which are depicted in 
figure 41: 


lOJ = M(s + ds) - M(s) + r^*(s) ds x F(s + ds) 


(s) ds 

+ X w(s) ds 


(A25) 


where 


0) = + co^e^ + 0)303 


M = + M2e2 + M3S3 


+Fe + Fe 
11 2^2 33 
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and I is the inertia tensor of the segment relative to the ( 02 ^ 02 ' ^ 3 ^ axis 
system. The principal axes of inertia are assumed to coincide with the 
( 61 , 62 / 03 ) axes so that the inertia tensor is diagonal in form. The diagonal 
elements of I are 


II - plj. 

ds 

(A26) 

^2 = P^p 

ds 

(A27) 

I 3 = PIz 

ds 

(A28) 


where Ip is the polar moment of the cross section and p is the mass density 
of the material of the ring when it is at rest. 

The components of the angular velocity 00 consist of the nominal spin 
rate projected onto the ( 6 - 1 ^, 62 / 03 ) triad and the angular rate of change 

of the triad with respect to the rotating ring: 


(jo-i = “^23 + 


R 

R + AR 


009 = w' +3 

^ R + AR 


03o = Q 


R 


r + Ar 


,(“■ - 1 ) 


(A29) 


(A30) 


(A31) 


The left side of equation (A25) may be written as 


= h + 1 


2 dt 


^2^2 . 


^3 dt^^3®3; 


(A32) 


Equations (A5) through (A7) and (A26) through (A31) are substituted into equa- 
tion (A32) , the indicated differentiations are performed, and products of the 
perturbation variables are dropped to yield a linearized equation for Ico in 
terms of the perturbation variables: 
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103 = 


(Er(5^ ”• - ») - (I. - i 

* ["p(® * i-^ “') " »r - "' - “ 


■ 5-TAi("' - sKjp '5= 


(A33) 


The first two terms of the right side of equation (A25) may be approxi- 
mated by 


M(s + ds) - M(s) 


9s 


Kn 


' 2“2 " '■' 3 '= 3 ) 


+ MoGo t M-^eo I ds 


(A34) 


Equations (A5) through (A7) are substituted into equation (A34) and the indi- 
cated differentiation is performed. Because M]^ and M 2 as given by equa- 
tions (A17) and (A24) consist of terms which are proportional to the perturba- 
tion variables w and 3, and M 3 as given by equation (A18) consists of 
terms which are proportional to the perturbation variables u and v, products 
of the perturbation variables and the quantities M]_, M 2 / ^1 * » ^2*^ ^3* 

are dropped for a linearized solution: 


M(s + ds) - M(s) 




+ 


1 

R + AR 






+ 



* i”i 



R 

R + AR 




ds 


(A35) 


Using the definition of 63 in equation (A5) and ignoring second order 
terms in the variable ds , the third term of the right side of equation (A25) 
may be written as 


rc* (s) ds X F(s + ds) 


kc* (S) I |F-:^ (s) 






Fi(s) 



ds 


(A36) 
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Equations (A5) through (A7) are substituted into equation (A36) and products of 
the perturbation variables are dropped to obtain 



(A37) 


Three equations which relate the internal forces and moments with inertia 
effects are obtained by substituting equations (A33 ) , (A35) / and (A37) into 

equation (A25) . Since the last term of the right side of equation (A25) is 
second order in the variable ds , it is dropped: 




R 


^\R + AR 


w - 


+ p(i^ - 


= Ml' - + 


7 — w* M, 

R + Ar / 3 


^ u\ . R_LM gj, 

V R R 3 R 1 


(A38) 




= ^2’ + 


S "1 ♦ (f - "")”3 ^ ”■'^1 - ("' - («3) 


-pi 


R 


2 R + Ar 


(=' -|)= V - ( 


S' - " s)s - 6S <«»> 


Substituting equations (Al7) , (A18) , and (A24) into equations (A38) 
through (A40) and using small-value approximations for the perturbation 
variables yield the following relations: 
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R + AR / R \ . 

r ^ ) w' 

R 3 rl R + AR / V 


3L\ + ( GJ + w 

R y I R + AR j^2 R + AR 


'r(s^“' - “S) " ® 


(A41) 


^ - V^VR + AR 


w* - 


GJR g„ _ J R . + ^ J \1. 

R + Ar ^Vr + Ar r ^ R z)^2 


+ T— GJ + E + 

R + Ar \ r 


-I fl 

R 


r + Ar 


(A42) 


R + Ar 




+ + V 

R\ R 




(A43) 


Equations (A41) and (A43) are formulas for the shear forces acting across the 
cross section as functions of the perturbation variables. Equation (A42) is 
the relationship between the angle of twist 3 and the out-of-plane 
deflection w. 


Equations of Motion 

Given a small segment of the ring, its motion obeys Newton's law which is 
described by 


r^pA ds = dF 


(A44) 


Substitution of equations (A4) and (A9) into equation (A44) yields a system of 
partial differential equations: 


lQi - 2Q-V - ( 


R + Ar + u)1 = W, + F, ' :^ + fu' - 

-I 1 1 R R + Ar V R/ ^ 


ir’ - Rf 1 ^ 
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Pa[v + 2«i - + Fj' + i - |)wi + WW 3 

+ ("" - t)''i - K'’’ - i )''2 + ""'■3] + I '■3 

pAw = W 3 + F^* + ^ ^”^ 2 ) “ ^^1 ” (A47) 


Solving the quiescent case of equation (A45) where the externally applied 
forces W and the perturbation variables are zero yields an expression for 
the equilibrium tension in the ring: 


F 


2 “ 



(A48) 


The equilibrium strain in the ring is obtained by equating equations (A48) 
and (A14) where the perturbation variables are set to zero: 


AR (Rfi)^ 



(A49) 


Let the externally applied force densities W-j^, W 2 / and be small 

perturbation quantities. Substitutions of equations (A14) , (A41), and (A43) 

into equations (A45 ) , (A46) , and (A47) and using equations (A42) and (A48) 

yield a system of partial differential equations of the perturbation variables 
for a steadily spinning ring with small external forces. In addition, if the 
assumptions 


e 


o 



1 


(Rfl)^ 




e 


O 


<< 



e 


o 



(A50) 

(A51) 

(A52) 

(A53) 
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are used and we note that Ip^ = then the resulting system of 

equations is 

W-, 

li - 2fiv - = _ _ o)g2g^2(j^4,j„.. ^ 2 r2u" + u) - 

+ f2^(R^u" - Rv' ) + S 2 ^(R^ii" - Rv') 

•• /-N • ^2 2 2 2 2 

V + 2f2u = (Ru‘ + R v") + Ru' + S (Rii' - v) 

p/\ s z 

W 3 

w = — - (r'^w"" - R^3") + (R^w" + r23") + q2j^2^„ 

+ Sr^R^(w" - 2J^3' + 

Sp^(3 + ^w') + - S^^^Q.(w' - fi3) = C0r^Sj2(R23" + Rw") 

+ Wg^S^2(Rw" - 3) 


(A54) 


(A55) 


(A56) 


(A57) 


where 



(A58) 

(A59) 

(AGO) 

(A61) 

(A62) 

(A63) 
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The terms 03^ and o)^ may be interpreted to be the speed of longitudinal 
and transverse sound waves in the material divided by the ring radius, and 
Sr/ Sp, and Sj are geometric characteristics which are analogous to the 
inverse of the slenderness ratio used in the design of columns (ref. 9) . 

Equations (A54) through (A57) form a system of partial differential equa- 
tions which describe the small-perturbation dynamic behavior of a steadily 
spinning slender ring which is influenced by small external forces. Equa- 
tions (A54) and (A55) describe the coupled radial and tangential motions for 
which the ring remains in its original plane. Equations (A56) and (A57) 
describe the coupled out-of-plane deflection and ring twisting motions. The 
out-of-plane and the in-plane motions are uncoupled from each other. 

Because of the absence of structural damping terms in the stress-strain 
relationships of equations (A12) and (A22) , the homogeneous solutions of the 
above equations are periodic in time. If damping were considered, the solu- 
tions would include exponentially decaying functions. It is expected that any 
practical ring would be very lightly damped so that the development of this 
paper would apply at least for short-term behavior. 


56 



APPENDIX B 


SOLUTION OF VIBRATION FREQUENCY EQUATION 

During the development of the characteristics of the structural dynamics 
of the ring performed in the main text, equations defining the vibration fre- 
quencies of the structural modes were found to be of the form 

f (0)) = OJ^ - ao)^ + bo) + c = 0 (Bl) 

where oa is the frequency of vibration; a, b, and c are positive real 
values; and b has a small value. If the roots of f ((jO) = 0 are real, then 
the motions of the ring are periodic in time. If any of the roots are complex, 
then the structural vibrations of the ring are unstable. This appendix pre- 
sents a procedure for calculating the roots of f {(jO) =0. A sufficient condi- 
tion for the existence of four real roots is also presented. 

Conditions on the values of a, b, and c in equation (Bl) are found so 
that if these conditions are satisfied, then there are no complex roots of 
f (tjo) = 0. Consider the special case of equation (Bl) where the value of b 
is zero: 

g (oo) - - ao)^ + c (B2) 

The roots of g (w) = 0 are easily found to be 



2 

These roots are all real and distinct if a > 4c . 

From equations (Bl) and (B2) , the roots of f (oo) = 0 satisfy the follow- 
ing equation : 


g (o)) = -bo) 


(B4) 


A typical curve for g (OJ) where a^ > 4c is presented in figure 42. The 
intersections of g (oo) and a straight line passing through the origin with 
slope -b occur at the roots of f (oo) = 0. For sufficiently small values of 
b, four intersections occur implying the existence of four distinct real roots. 
For sufficiently large values of b, only two intersections occur implying the 
existence of two real roots and two complex roots. The line *“h^^oo in fig- 
ure 42 has two intersections and is tangent at a third point on the g (oJ) 
curve. In this case, there exists a double real root at the tangent. For 
b > b^, there exist complex roots, and for b < b^, there exist only real roots 
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of f (co) = 0. Therefore, the value of is the greatest lower bound on the 

values of b for which there are complex roots of f ((jO) = 0. 

Consider the line -bj^oo in figure 42 which passes through an extremum 
of g(oo). The value of b^ is given by 


b 


L 


- 4c 

\j 8a 


(B5) 


The value of bj^ being less than b^ is a lower bound on the values of b 
for which there may be complex roots of f ((jo) = 0. Therefore, if 


a^ > 4c 


and 


b S 


a^ - 4c 


(B6) 


(B7) 


2 

then all roots are real and distinct. If b is nonzero and a = 4c or if 
a^ < 4c, then there exist complex roots of equation (Bl) . If equation (B6) is 
satisfied but equation (B7) is not, then a more detailed analysis is required 
to determine whether the roots are all real. 

The roots of the quartic polynomial equation f (co) = 0 can be calculated 
by algebraic techniques, for example, Ferrari ' s method (ref. 16). However, a 
different procedure is used in the following analysis. 

Given that f (cjo) has four real roots, then it can be factored into the 

form 


f (oo) = (0) - (jOq “ COjl) (00 - OJq + 00j_) (00 + OOq - OO 2 ) (00 + OOq + OO 2 ) (B8) 

where oOq, 00 ^, and 0 O 2 are real numbers. 

Multiplying out equation (B8) and equating coefficients with equation (Bl) 
yield a set of simultaneous nonlinear equations in oOq, oo^^, and 0 O 2 : 


00 


+ 00. 


+ 200 


0 


(B9) 


2(002^ - 00 ^ 2 ^ OOq = b 

(“0" - - "2 




= C 


(BIO) 

(BID 
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2 

Solving equation (B9) for 0 O 2 and substituting the result into equation (BIO) 
yield an equation for as a function of 03 q: 


0)1 


2 





(B12) 


2 

Similarly, solving equation (B9) for o)^^ 

2 

equation (BIO) yield an equation for 0)2 


and substituting the result into 
as a function of 0 )q: 


0) 


2 


2 




(B13) 


Substitution of equations (B12) and (B13) into equation (Bll) yields a poly- 
nomial in 0 )q! 


0)o 


6 



(B14) 


The roots of this equation may be found by application of Tartaglia*s method 
(ref. 16) . In this paper only those roots of equation (B14) with the smallest 
magnitude are considered. Assuming that the magnitudes of these roots are much 
smaller than the magnitudes of the remaining roots, an approximate solution may 
be obtained by ignoring the first term of equation (B14) : 


■ (B15) 

(a2 - 4c) + \J (a^ - 4c) ^ - Sab^ 

For more precision, an iterative procedure may be used for solving equa- 
tion (B14) by substituting an estimate (or previous solution) for for 

the first term of the equation and solving as if it were part of the constant 
term as follows: 



(a^ - 4c) + 


Ip - 

- 4c) ^ - Sajb^ - 


(B16) 


where the subscript i denotes the ith solution of the iteration and 

(“o^)o - “■ 

without loss of generality, only positive values of o)q are considered. 
With this condition, examination of equations (Bl2) and (B13) shows that the 
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magnitude of ( 1)2 is larger than the magnitude of Solving equation (B9) 

for 0)2^ and substituting the result into equation (Bll) yield a polynomial 

in 

+ To + WQ^^a - 303q2^ = 0 (B17) 


Similarly solving equation (B9) for substituting the result into 

equation (Bll) yield an identical polynomial in 002^’ Equations (B12) and (B13) 
show that there are only two solutions for 03^ and two solutions for 032 . 

Since 1^2 1 ^ 1^1 I' larger roots of equation (B17) are values of 032 

and the two smaller roots are values of 03j_ : 



(B19) 


From equation (B8) ^ the four roots of f (o3) = 0 where f (03) is given by 

equation (Bl) are 03 = 03q ± 03^^, -“03q ± 0)2 where 03q is calculated by the 

approximate equation (Bl4) or the iteration equation (B16) , 0)2 is calculated 

by equation (Bl8) , and 03j_ is calculated by equation (B19) . All the roots are 

real if the conditions of equations (B6) and (B7) are satisfied. However, 

there are cases which have all real roots but which violate the condition of 
equation (B7) . 
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TABLE I 


CHARACTERISTICS OF EXAMPLE RING 


Symbol 

Value 

Symbol 

Value 

R 

360 m 


1.12 X 10“^ 

A 

1.96 cm^ 

Sr 

in 

1 

o 

t — 1 

X 

OJ 
1 — 1 

1 — i 


0.320 

Sp 

1.59 X 10“5 


0.320 

Sj 

1.46 X 10"5 


0,640 cm"^ 


1.52 

J 

0.538 cm^ 


12.5 km/sec 

P 

1 . 8 g/cm^ 


7.8 km/sec 

E 

280 GPa 


34.7 rad/sec 

G 

110 GPa 


21.7 rad/sec 











Spin direction 




Figure 10.- Out-of-plane standing wave motion for 
high spin rate. 










Figure 17,- In-plane progressive bending mode for low spin rate, 

Q « j = 2. 



‘igure 18.- In-plane regressive bending mode for low spin rate. 
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Figure 19.- In-plane bending standing wave as seen by 
a fixed observer for low spin rate. 0 . « S^O) . 
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Figure 23 In-plane bending standing wave as seen by a fixed 
observer for high spin rate. << << 03g. 




Figure 25.- In-plane compression standing wave 
as seen by fixed observer. 
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Figure 28.- Out-of-plane modal frequencies. 
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Figure 31.- Out-of-plane modal velocities versus spin rate 
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Figure 33.- In-plane modal frequencies. 
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Figure 34.- Ratio of tangential motion to radial motion. Q = 1 rad/sec 
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Figure 35.- In-plane nodal velocities. 






Figure 36.- In-plane modal velocities versus mode number 
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Figure 40.- Curvilinear coordinates. 
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Figure 41.- Short ring segment with forces and 
moments acting upon it. 
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